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Abstract 



In this work, we calculate the branching ratios and CP asymmetries of 
the decays of B — > tttt and ttK with the frame of QCD factorization in the 
heavy quark limit. We also compare the results with the estimates by using 
generalized factorization and experimental measurements. 
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In the past two years, the CLEO collaborations had improved their measurements for 
the branching ratios of B — > tttt and ttK for several times. The latest results of branching 
ratios of these modes are: 
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These two-body non-leptonic charmless B decay modes play a very important role in studying 
CP violation and the heavy flavor physics. Theoretically, due to our ignorance of hadroniza- 
tion, it is difficult to relate CP violations with the parameters of fundamental theory. Ef- 
fective Hamiltonian is our basic tool, but we do not know how to calculate hadronic matrix 
element from first principles. Conventionally we resort to naive factorization assumption [|TJ , 
which expresses the hadronic matrix element in terms of various meson decay constants and 
form factors. However, it is well known that non-factorizable contribution really exists and 
can not be neglected numerically. To remedy factorization hypothesis, Ref. |5]|| introduced 
a phenomeno logical parameter N e ff, which is commonly called generalized factorization. 
Because in principle N e ff is process-dependent, it is still not a satisfactory approach. 

In last year, Beneke et al. gave a NLO calculation of B — > nix in heavy quark limit. 
In this limit, the hadronic matrix elements for B — >• ttti can expanded by the powers of a s 
and AQ C D/ m b as follows: 

(nn\Q\B) = (7r|ji||S)<7r|j 2 |0) • [1 + £r n < + 0(A QCD /m b )}, (2) 

where Q is a local four quark operator in the weak effective Hamiltonian and j\,2 are bilinear 
quark currents. Neglecting the power contribution of they pointed out that in 

the heavy quark limit the radiative corrections at the order of a s can be calculable with 
PQCD method. Furthermore, the 'non-factorizable' contributions from hard scattering with 
spectator quark in B meson can also be calculable within the frame of PQCD. Then the 
decay amplitude can be expressed by the convolutions of the hard-scattering kernels and 
several light-cone distribution amplitudes of the mesons. So all of these can be summarized 
into a factorization formula for B — >• tttt as follows: 

l l 

(7T(p')TT(q)\Q t \B(p)) = F B ^(q 2 ) J dxT t ? (x)$ n (x) + J d^dxdyT H ^,x,y)^ B (0%(x)%(y). 



(3) 

We call this factorization formalism as QCD factorization. In the above formula, $b(£) 
and $7r(a;) are the leading- twist light-cone distribution amplitudes of B and pion mesons 
respectly, and the T?' 11 denote hard-scattering kernels which are calculable in perturbative 
theory. Beneke et al. regarded the spectator quark as a soft quark translated to one 
pion in final state unless it undergoes a hard interaction, so the transition form factor 
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F B ^ n (q 2 ) can not be calculated in the perturbative frame, and it will be survived as a 
free nonperturbative parameter in QCD factorization. T/ at the order of a s includes two 
topological classes, one is hard gluon scattering like vertex corrections which is described by 
the first four diagrams in Fig.l, and we call it as hard-scattering corrections below; the other 
is called penguin correction, which is described by diag.(e) and (f) in Fig.l. T* denotes 
the hard spectator scattering contribution which is described by the last two diagrams in 
Fig.l. In this frame, they neglected W-annihilation and W-exchange topologies, so they 
show that 'non-factorizable' contributions in the conventional factorization are calculable in 
heavy quark limit. Then we do not need to employ a phenomenological parameter N e ff to 
compensate the 'non-factorizable' effects in principle. Thus we will apply this approach to 
calculate the branching ratios and CP asymmetries for the decays B — > -kit and 7rK in this 
paper. Furthermore, since the effects of electroweak penguins can not be neglected in some 
decay modes of B — > ttK, we will use full effective weak Hamiltonian including electroweak 
penguin operators and add some one-loop QED penguin correction to the calculation of 
hadronic matrix elements of effective operators which can also be described by diag. (e) and 
diag.(f) (replacing gluon with photon) in Fig.l in our computation. 
The \AB\ = 1 effective Hamiltonian is given by JK| 



Heff ~ — 



V2 
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E v q CMQKri + cmqpm + £ cmqm 

J=u,c \ k=3 J 



+ h.c, (4) 



where v g = V^V^for b — > d transition) or v q = V q t,V* s (for b — > s transition) and Cj(//) are 
Wilson coefficients which have been evaluated to next-to-leading order approximation. In 
the Eq.(l), the four-quark operators Qi are given by 

Ql = {u a bp) v ^ A (qpu a ) V - A Qi = (c a bp) v ^A(qpc a ) V -A 

Q2 = (u a b a )v-A(q~i3Up)v-A Q\ = (c a b a ) V -A(qpcp)v-A 

Qi = (q a ba)v-AY,{qpq'p)v-A Qi = {qpb a ) V -AY.{q a qf})v-A 

q' q' 

Q5 = (qab a )v~AY,(qaqi3)v+A Q& = {qpb a )v-AY,{q a qR)v+A fK \ 

ql q ' ^ ' 

Qi = l{qab a )v-AJ2^q'(qi3qi3)v+A Qs = l(q/3b a )v-AJ2 e q'{.q a qp)v+A 

q' q' 

Q9 = \{qab a )v~AY, e q'{qi3qp)v-A Qw = \{qpb a )v-AY, e q'{q a qp)v-A 

q' " q' 



with Q\ and Ql being the tree operators, Q 3 — Q 6 the QCD penguin operators and Q 7 — Q w 
the electroweak penguin operators. With the renormalization group method, we can evolve 
the Wilson coefficients Cj(/i) from the scale \i = mw to /1 ~ tub- Because we will give a 
NLO calculation here, b — > s(d)g and b — > s(c?)7 effective Hamiltonian must be included. 
They are 

Q 

1-L'eff = --J=V t \ C 7~iQ7l + CsgQsg] + h.C, (6) 

where 

Qi~f = ^ m bq a ^ iu {l + ^ 5 )b a F^, Qsg = ^-m b q a a^t a al3 bfjGl u , (q = d or s). (7) 
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Following the method of Ref. 0, we try to evaluate the matrix elements of Qi to the 
order of a s and ot em . In fact, this work is to calculate hadronic matrix elements of Qi to one 
loop. In quark level, Ref |3J] gives an expression 

(Q) = [i + ^m s + ^m e ] ■ (Q) . (8) 

47T 47T 

Here rh s and m e represent the one loop corrections of QCD and QED respectly. These cor- 
rections are divided into two classes. One is hard-scattering correction, the other is penguin 
correction. While the external quarks are on mass-shells, the hard-scattering corrections 
bring infrared divergences which can not vanish after summing over all this kind of pertur- 
bative diagrams no matter with gluon or photon exchange. However, in hadron level, with 
QCD factorization in the heavy quark limit, these infrared divergences from hard-scattering 
exchanging with gluon can be canceled after summing over all this kind of perturbative di- 
agrams (the first four diagrams in Fig.l). This has been shown in Ref. ||. But for the case 
of QED hard-scattering (exchanging photon), the infrared divergences can not be canceled 
after summing over all diagrams even in the heavy quark limit. It is because that the color 
structure of QED hard-scattering corrections is different from that of QCD. We expect that 
those divergences could be canceled by soft photon radiative corrections. However, nobody 
knows how to include these radiative corrections in exclusive decay channels. As a em is very 
small, we will neglect QED hard-scattering corrections in our computation. On the other 
hand, the penguin corrections are calculable not only for the case of QCD but also for that 
of QED. In quark level, considering only the contributions from tree operators, the penguin 
corrections can be written in NDR scheme as ||] 

1 2 

(m s (/i))i3 = (m s (/i))i5 = — [-- + G(m q , q, //)], (9) 

1 2 

(rh s (ii))u = (m,(^))i 6 = -^hg + G ( m qi<l, v)], ( 10 ) 

4 2 

(m e ( / u)) 27 = (m e (/i)) 2 g = -ghg + G ( m <?> 9, ft)], ( n ) 

4 2 

(m e (/i)) 17 = (m e (/i))i 9 = - g[-g +G{m q ,q,ti)], (12) 

with 

} m 2 — x(\ — x)q 2 — ie 
G(m q , q, /j,) — — 4 I dx x(l — x) In — . (13) 



o 



Here, q 2 in G(m q , q, /i) remains uncertain for the calculation of exclusive B decays. How- 

2 

ever, if we take q 2 around in computations with conventional factorization scheme, this 
variation does not change the results too much 0. With the frame of QCD factorization 
in the heavy quark limit, there will be no uncertainty for q 2 , this will be shown below. So, 
similar to what done in Ref. |2j, we will take both QCD and QED penguin corrections into 
account. Surely, for the importance of electroweak penguins in the decays B — > ttK, we 
must take QED penguin corrections into account as well. 

Then in heavy quark limit, the amplitude for the decay of B to two light pseudoscalar 
mesons P and P' can be written as follows: 
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A(B -> PP>) = % E E V^Pmi^F, (14) 

V ^ p=u,c i=l,10 

where v p is CKM factor and (PP'\Qi\B) F is the factorized matrix elements which can be 
expressed by the product of the relevant decay constant and form factor. Taking a\ = a c 2 = 
and assuming the distribution amplitudes of light pseudoscalar mesons symmetric, we obtain 
the QCD coefficients a\ at next-to-leading order (NLO) in naive dimension regularization 
(NDR) scheme (except and a 8 which read at leading order here for some special reasons). 
In order to cancel the scheme dependence in our calculation, we must take Wilson coefficients 
Ci in NDR scheme as well. Then the explicit formulas of a\ can be written as 

(15) 

- = ^ + t + S^ < 17 > 

2 SO 2 f 

+ (Gp(s p ) - -)d + (G P (0) + G P (1) -f P - fj, 1 + y )C 3 - -j-C, 

+(3G P (0) + G P ( Sc ) + G P (1))(C 4 + C 6 ) + Gp i8 C 8G ], (18) 
a 5 = C 5 + | + ^fc 6 (-F-12), (19) 

a 6 = C 6 + ^, (20) 

«r = C 7 + ^ + ^C 8 (-F-12), (21) 

as = C 8 + ^, (22) 

a^C 9 + f + ^C 10 F, (23) 

+(§ " G P (s p ))(Ci + ^) + (\- Gp(s c ) + ^)(C 3 + §) 

+(-Gp(0 + ^)(C 5 + ^) - ^C 77 Gp, 8 ]. (24) 

Here iV = 3 (/ = 5) is the number of color (flavor), and Cp = is the factor of color. 

We define the symbols in the above expressions as the same as Beneke's, which are 

F = -12 In -^-18 + + (25) 
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f P = J dx g(x)$ P (x), Gp, 8 = J dxG s (x)$p(x), (26) 



1 

G P (s)= f dx G(s,x)$ P (x), (27) 



with 



1 — 2x 2 
g(x) = 3— lnx — 3iii, G 8 (x) = , (28) 
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G{s,x) = -aJ duu{l-u)\n{s-u{l-u){l-x)-ie). (29) 
o 

Here <&p(x) is the leading twist distribution amplitude of the light meson, and s$ = mf/mf . 
(rrii is the mass of quark appearing in the penguin loop.) The contribution from the hard 
spectator scattering are reduced into the factor fp 1 which is written as 



In above expression, fp> (fp) is the pseudoscalar meson (B meson) decay constant, mp the 

B meson mass, F+~* p (0) the B — > P' form factor at zero momentum transfer, and £ the 

light-cone momentum fraction of the spectator in the B meson. 

One can find that our expressions are a little bit different from those in Ref . @ . We think 

that there is an extra term — ^C^ in Eq.(18) which is missed in Eq.(8) of Ref. [Q], but this 

difference is not large enough to make the numerical results change too much. In addition, 

in the modes B — > irK, fp 1 should be changed with the corresponding meson which contains 

the spectator quark. But numerically the ratio gfes is nearly equal to Trfer, so fU ~ fl 1 , 

p + * + 

we will not distinguish them in computation below. 

For the coefficients a 6 and as, we want to give some comments. One can see that the 
singularity of the hard-scattering function Gs{x) is at the endpoint x — 1. If the distribution 
amplitudes did not fall off fast enough at the endpoints to suppress these singularities in hard- 
scattering functions, our calculations would not be infrared safe. For the case of a!_ 5 and 
a 7,9,i0; we might use the twist-2 distribution amplitudes for pion or kaon which can cancel 
the infrared divergences. Unfortunately, for the case of a 6j8 , we must employ the twist-3 
distribution amplitudes. As pointed in Ref. p|,|l2|, the twist-3 distribution amplitudes for 



pion and kaon do not fall off at endpoints. So if this was true, the infrared divergences could 
not cancel in the calculation. Therefore, as mentioned at the end of Ref. , the factorization 
formula breaks down in this case. We also noticed that in recent papers ||||], the authors 
employed the twist-3 distribution amplitudes for pion and kaon which fall off fast enough 
at the endpoints, but that is just a model and not a prediction from QCD sum rule or 
other non-perturbative approaches. So we still employ the asymptotic twist-3 distribution 
amplitudes of pion and kaon as $p(x) = 1, then the singularities at the endpoints in hard- 
scattering kernels can not be suppressed by the distribution amplitudes in the case of a 6 and 
as- Like what Beneke et al. do, we take a§ and a§ at leading order here. 
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In the B rest frame, the two body decay width is 



where 



T(B^PP') = — \A(B ^ PP')\ 2 ^l, (31) 
8% m B 



= [(m 2 B - (m P + m pi f)(m 2 B - (m P - m P ,) 2 )} 2 

1m b 

is the magnitude of the momentum of the particle P or P' . The corresponding branching 
ratio is given by 

Y(B -> PP') 

BR(B -> PP 1 ) = — I (33) 

The direct CP asymmetry Acp for B meson decays into PP' is defined as 

= |A(j?^ppQ| 2 -|^^pp)| 2 

CP \A(B^PP')\ 2 + \A(B^PP')\ 2 ' 1 J 

Here we do not consider B° — B° mixing just for simplification. 

Because the momentum fraction distribution of the spectator quark in B meson is peaked 
at Aqcd/'itT'B, we will take the distribution amplitude for B meson as 

MO = (35) 



^ is equal to the parameter m B j\ B m Ref. ||. Here we take e B = 0.05. For simplification, 
we will take the asymptotic form for the leading twist distribution amplitudes of pion and 
kaon as same as Ref ||: 

^{x) =6x{l-x), $ K (v) = Gv(l-v). (36) 

After straightforward calculations, we carry out the branching ratios and direct CP 
asymmetries for the decays B — > mi and ttK at two different renormalization scales fi = 
5.0 GeV and 2.5 GeV which are listed in Tab. 2 and Tab. 3 respectly. (Wilson coefficients 
in NDR scheme [|6|,[7l,|10] at two scales are listed in Tab.l.) As a comparison, we also show 
the results with the conventional factorization (BSW factorization and N e ff = 2) in the 
last two tables. In our computation, we take the Wolfenstein parameters of CKM matrix 
as follows: A = 0.8, A = 0.22, p = —0.12, r\ = 0.34. The corresponding decay constants 
and form factors are taken as follows §: f n = 0.13 GeV, f K = 0.16 GeV; F B ^ 7r (0) = 0.33 
and F B ^ K (0) = 0.38. Quark masses are taken as: m^ = 4.8 GeV, m c = 1.4 GeV and 
m s = 105 MeV. 

From both Tab. 2 and Tab. 3, we find that QCD factorization enhance the contributions 
from penguins so much that the branching ratios of B — » nK become larger nearly by a 
factor of 2 than the estimate of the generalized factorization. It is because that the QCD 
and QED corrections are constructive to the amplitudes with the generalized factorization 
and they enhance the branching ratios. One also see that the contributions from electroweak 
penguins can not be neglected in the decays B~ — > ir°K~ and B® — > tt°K no matter with the 
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generalized factorization or QCD factorization, because the amplitudes of these two modes 
contain a term of the coefficient ag which is large enough to compare with the coefficients of 
QCD penguins. For the case of B® — > 7r + 7r _ , the differences between these two factorization 
schemes are not too apparent because that this decay mode is dominated by the coefficient 
a\. The radiative correction to a\ at the order of a s is very small comparing with the 
leading order part of a\. This is similar to that the effective coefficient a\ = C\ + C2/N e ff 
in generalized factorization is unsensitive to the variation of the phenomenological parameter 
N e ff. But for the coefficient 02, it is very different for its small leading order part. So for the 
modes B~ — > 7r°7r~ and — > 7r°7r°, the differences between these two factorization schemes 
are very large. 

For the case of Acp, the differences between generalized factorization and QCD factor- 
ization in the heavy quark limit are quite large because in QCD factorization the imaginary 
parts enter the decay amplitudes through hard scattering kernels contribution. Then the 
strong phases of some aj (i=even number) dominant processes may be changed very much, 
then CP asymmetries of these modes can be dramatically large. One can see that CP asym- 
metries of B® —>■ 7r°7r° is about 80%, but this magnificent direct CP asymmetry is very hard 
to observe for its small branching ratio. For the modes B — > ttK , the CP asymmetries change 
much as well. In the modes of B — > ttK, the signs and magnitudes of CP asymmetries are 
changed comparing with the results of the generalized factorization. 

The authors of Ref. @ pointed out that the amplitudes derived from the QCD factor- 
ization in the heavy quark limit are independent of the renormalization scale physically. 
Numerically, we still find that the dependences in our results of the branching ratios for 
the decays B — > ttK are visible. Comparing the results at the scale fi = 5.0 GeV and 
/i = 2.5 GeV, this variation brings about ±20% uncertainty to the estimates for the branch- 
ing ratios of B — > irK. As shown in Ref. and our paper, the scale dependences of the 
results B — > nn are small. In recent calculation of B — > D^tt" in the heavy quark limit 
]T3| , the scale dependences are also small. Maybe for the case of pure tree or tree dominant 
processes, the computation with the frame of the QCD factorization cancels the dependence 
of the renormalization scale very well. But for some pure penguin or penguin dominant 
processes, the scale dependences are visible. However, these dependences are smaller than 
the estimates with the generalized factorization. 

We also show the dependences of the branching ratios and direct CP asymmetries on the 
weak phase 7 = argV^ respectly in Fig. 2 and Fig.3. In both Fig. 2 and Fig.3, the results of 
(a) and (b) are carried out with the QCD factorization in the heavy quark limit at the scale 
H = 2.5 GeV. 

From Fig.2 and Fig3, we find that the results are in favor of the experimental measure- 
ments when 90° < 7 < 270° because the branching ratio of B® — > 7r + 7r _ becomes small 
and the branching ratios of B~ — > tt°K~ , B~ — ► tt~K° and B® — > n + K~ are closer in that 
region. This is consistent with the fit for 7 of CLEO and other researchers PJT4} . 

We note that our estimate of branching ratio of B® — > 7r + 7r _ seems larger than the 
experimental measurement even if we take 7 > 90°. And the branching ratio of B® — > tt°K° 
is about 3 times smaller than the central value of experimental result and unsensitive to the 
variation of 7. Due to uncertainties of the form factors and F+~* K , we try to vary 

these form factors in a relevant region to make the branching ratio B^ — > 7r + 7r~ smaller while 
that of Bj — > tt°K° is larger. We find the branching ratios of B — ► tttt and ttK are rather 



8 



sensitive to the form factor F+^ n than F^ K . And both of them increase with F^ w . So 
our attempt faces a failure. But now the errors in present measurements of CLEO are still 
large, and some uncertainties remain in theoretical frame, such as, the light-cone distribution 
amplitudes of the mesons, the heavy to light transition form factors and etc. Therefore, the 
disagreement between our prediction and the present experimental measurement is not so 
significant. It needs us more detailed study with the improved experimental measurements 
and theoretical approaches in future. 

In conclusion, QCD factorization can give an estimate of strong phases in B charmless 
decays from final state hard scattering. This might be beneficial to extracting the weak 
phases from CP asymmetries in B decays. But due to the theoretical uncertainties such 
as meson light-cone distribution amplitudes and heavy to light transition form factors, the 
prediction for branching ratios within the frame of QCD factorization remains a little bit 
ambiguous. On the other hand, since that is not a very large scale, maybe the corrections 
at the order of Agc^/m^ are needed. So, how to develop a complete perturbative frame of 
heavy quark expansion in heavy to light decays will be a potentially interesting and beneficial 
work. 
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TABLES 





Ci c 2 c 3 c 4 c 5 c 6 


H = 5.0 GeV 


1.150 -0.313 0.016 -0.033 0.009 -0.042 


H = 2.5 GeV 


1.117 -0.257 0.017 -0.044 0.011 -0.056 






CV C 8 C 9 Cio cv 7 C 8G 


H = 5.0 GeV 


-2 x 10~ 5 38 x 10" 5 -0.0103 0.0021 -0.300 -0.144 


H = 2.5 GeV 


-1 x 10" 5 50 x 10~ 5 -0.010 0.002 -0.336 -0.158 


TABLE I. Wilson coefficients in NDR scheme 


H = 5.0 GeV 
Decay Modes 


Branching Ratios CP Asymmetries 
Generalized FA QCD FA Generalized FA QCD FA 


B~ - 7T°7T- 

B° d -► tt+vt- 

5° -► 7T ^ 


6.33(6.71) 4.30(4.60) -0.1(0) 0(0) 
7.36(7.45) 7.93(7.97) 2.4(2.5) -3.2(-3.6) 
0.38(0.49) 0.14(0.21) -7.5(-5.7) 78.4(64.7) 


B" - 7T°K- 
B" -»• 7T-i?° 
BO _ 

B° _ ^K 


6.00(3.81) 10.2(7.10) -4.5(-7.1) 4.1(5.7) 
5.50(5.81) 12.3(11.1) -O.l(-O.l) 0.4(0.4) 
7.80(7.01) 15.0(13.9) -5.4(-6.1) 2.4(3.0) 
1.47(2.7) 3.90(5.20) 3.8(4.0) -3.2(-2.8) 



TABLE II. In this table, branching ratios are in the unit of 10 -6 and CP asymmetries are in 
the unit of one percent. And the values in the brackets are the results without considering the 
contributions of EW penguins. 
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fl — Z.O IjreV 




Branching 


Ratios 


CP Asymmet 


nes 


Decay Modes 


Generalized FA 


yL/JJ rA 


Generalized FA 


yL/JJ rA 


D— „0„ — 

B u — > 7T 7T 


6.07(7.08) 


A At ( A 

4.41(4.65) 


-0.1(0) 


0(0) 


5° tt+vt- 


6.94(7.02) 


7.54(7.50) 


2.8(2.8) 


-3.2(-3.8) 


5° - 7r°7r° 


0.63(0.77) 


0.26(0.37) 


-5.4(-4.4) 


73.4(61.1) 




8.50(5.89) 


11.9(9.30) 


-3.2(-4.7) 


4.3(5.3) 


B- -»• vr^i? 


9.08(9.44) 


15.2(15.1) 


-O.l(-O.l) 


0.4(0.4) 


£° - 7T+K- 


11.90(11.07) 


17.9(18.1) 


-3.5(-3.9) 


1.9(2.2) 


B° d - 7T ^ 


2.78(4.41) 


5.00(7.00) 


2.2(1.3) 


-3.8(-3.2) 



TABLE III. In this table, branching ratios are in the unit of 10 -6 and CP asymmetries are in 
the unit of one percent. And the values in the brackets are the results without considering the 
contributions of EW penguins. 
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FIGURES 




(e) (f) (g) (h) 



FIG. 1. Order of a s corrections to hard-scattering kernels T 1 and T n . The upward quark lines 
represent the ejected quark pairs from b quark weak decays. 




FIG. 2. Br(B tttt) x 10 6 vs. 7 = arg V* b (Figl.(a)) and Br(B -► ttK) x 10 6 vs. 7 (Fig.(b)). 
In Fig. (a), solid, dashed curves are for B® — > 7r+7r~ and 5,7 — ► 7r°7r~ respectly; in Fig.(b), solid, 
dashed, dotted and dot-dashed curves are for B~ — ► ir°K~, ir~K° and — > 7r + K~, 7r°K° 
respectly. 
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FIG. 3. A C p(B -► tttt) x 10 2 vs. 7 = aigV* b (Figl.(a)) and A CP {B -► vriT) x 10 2 vs. 7 
(Fig.(b)). In Fig. (a), solid, dashed curves are for B® — > 7r + 7r~ and — > 7r°7r _ respectly; in 
Fig.(b), solid, dashed, dotted and dot-dashed curves are for I?~ — ► ir°K~, ir~K° and 5° — > 7r + i"C~, 
7r°E"° respectly. 
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